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DESCRIPTION 

CHANNEL SIMULATION METHOD AND CHANNEL SIMULATOR 

Technical Field 
5 The present invention relates to a channel 

simulation method to simulate wireless channels in 
mul t i - ant enna communications, and to a channel simulator 
to simulate wireless channels in multi-antenna 
communications to promote development of wireless 
10 apparatus . 



Background Art 

Conventionally, in developing a cellular telephone, 
a base station for the telephone, and MT and AP of wireless 

15 LAN , there is needed an apparatus for simulating a wireless 
channel as an environment for performance evaluation of 
a development apparatus, i.e. a channel simulator. 

By comparing transmission characteristics obtained 
by providing signals transmitted from a development 

20 apparatus with simulated fading and receiver noise using 
a channel simulator with theoretical values or computer 
simulation values, it is possible to judge whether the 
development apparatus performs desired operation. 
Further, by reproducing a channel state of a running 

25 experimental course, it is possible to analyze a failure 
of a development apparatus occurring under an actual 
propagation environment. By thus using a channel 



, flCterlstlC evaluations of a development 

simulator, cnara 

apparatus can be performed in a room with ease. 

PIG l shows an example o£ a configuration of a 
conventional channel simulator. In channel simulator 
, 10 a transmission si g na! output from a transmission 
sy ; te m of a development apparatus ,0 is passed through 
a multipath channel configure* according to set 
parameters from control apparatus 30. At this point, 
a signal passed through each path is given the amplitude 
0 variation and phase variation .hereinafter, referred to 
, s onannel variations, simulating fading, and weighted 
uith a gain of each path to he added. The signal given 
the channel variations and receiver noise by channel 
simulator 10 is received and demodulated in reception 
„ system 50 of the development apparatus, and the 
demodulated signal is output to error rate measurement 

, vn Thus, it is possible to evaluate 

instrument 70. Tnus, 

performance of transmission system 40 and reception 
system 50 of the development apparatus by observing error 
20 rate measurement results obtained by adding various 
channel variations using channel simulator 10. 

A specific configuration of channel simulator 

■h.rt below Channel simulator 10 is 
will be described below. 

■ i„„ svstem 40 of the development 
connected to transmission system 

25 apparatus comprised of digital baseband processing 
section .digital BB processing section, 41, analog 
haseband processing section .analog BB processing 



/ 
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section) 42 and radio circuit 43, while being connected 
to reception system 50 of the development apparatus 
comprised of radio circuit 53, analog BB processing 
section 52 and digital BB processing section 51. In 
5 addition, in FIG.l, each line represents two baseband 
signal lines comprised of an I channel (in-phase i.e. 
the real part of a complex number) and a Q channel 
(quadrature i.e. the imaginary part of the complex number ) , 
except lines for connecting between radio circuits 43 

10 and 11 and between radio circuits 20 and 53. 

Digital data generated in data generator 60 is input 
to channel simulator 10 via digital BB processing section 
41, analog BB processing section 42 and radio circuit 
43 of transmission system 40. Digital BB processing 

15 section 41 is a portion to perform digital modulation, 
spreading and other processing when transmission system 
40 of the development apparatus is a C DMA (Code Division 
Multiple Access) transmission apparatus, while being a 
portion to perform digital modulation, inverse Fourier 

20 transform and other processing when the system 40 is an 
OFDM transmission apparatus. Analog BB processing 
section 42 is a digital/analog conversion circuit, and 
radio circuit 43 is a portion to perform upcon ve r t ing , 
signal amplification, etc. 

25 Channel simulator 10 has radio circuit 11 that 

performs processing inverse to that in radio circuit 43, 
i.e. processing such as downconvert ing , and analog BB 
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processing section 12 comprised of an analog/digital 
conversion circuit, and converts a signal from 
transmission system 40 into a digital baseband signal 
using radio circuit 11 and analog BB processing section 
5 12 . 

The digital baseband signal is input to multipath 
signal generating section 13 comprised of shift register 
14 and selector 15, and becomes a multipath signal in 
multipath signal generating section 13. More 

10 specifically, shift register 14 shifts the digital 
baseband signal input thereto by time obtained by dividing 
the maximum delay time of a path by a sampling cycle of 
analog BB processing section 12. 

Selector 15 selects signals corresponding to the 

15 number of paths from among signals output from each shift 
stage of shift register 14. At this point, multipath 
instruction signal SI indicative of the number of paths 
and delay time of each path designated from control 
apparatus 30 is input to multipath generating section 

20 13, and shift register 14 and selector 15 operate based 
on multipath instruction signal SI. In this way, 
selector 15 in multipath generating section 13 outputs 
a signal corresponding to each path under multipath 
environments . 

25 The signal corresponding to each path is output to 

a respective one of complex multipliers Al to Ak in 
ins t ant aneous variation ( Ra y 1 e i gh f ad i ng ) adding s e c t i on 
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16. Complex multipliers Al to Ak are supplied with 
complex gaussiannoises generatedby band-1 imited complex 
gaussian noise generating sections (LGN) Dl to Dk, 
respectively. In addition, each of band- limited complex 
5 gaussian noise generating sections (LGN) Dl to Dk is 
comprised of a white gaussian noise generating section 
and a Doppler filter, and generates a white gaussian noise 
limited in band to a range of maximum Doppler frequency 
S2 input from control apparatus 30. By this means, 

10 complex multipliers Al to Ak output respective signals 
of paths each provided with the instantaneous variation. 

The signals of paths provided with instantaneous 
variations are output to a plurality of complex 
multipliers, Bl to Bk, forming short-term variation 

15 adding section 17. Each of complex multipliers Bl to 
Bk is supplied with complex gain S3 corresponding to the 
respective path designated from control apparatus 30, 
and thus, short-term variation adding section 17 outputs 
a signal of each path provided with shadowing and distance 

20 variation. By this means, in channel simulator 10, for 
each path, a signal is formed which is provided with the 
instantaneous variation, shadowing and distance 
variation designated from control apparatus 30, and the 
signal of each path is all added in adders CI, C2,..., 

25 wherebya mul t ipa t h signal is fo rmed in which are reflected 
channel variations . 

The multipath signal is supplied to adder C3. Adder 
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C3 is also supplied with the white gaussian noise that 
is generated in white gaussian noise generating section 
(WGN) 21 and amplified in amplifier 22 to noise level 
S4 designated from control apparatus 30. In this way, 
5 adder C3 adds a receiver noise to the multipath signal. 

Analog BB processing section 19 and radio circuit 
20 respectively have the same configurations as those 
of analog BB processing section 42 and radio circuit 43 
of transmission system 40, perform digital / analog 

10 conversion on a digital BB signal given the channel 
variations and receiver noise, and then, further perform 
radio processing such as upconverting and amplification 
on the converted signal. 

An output signal of channel simulator 10 is input 

15 to radio circuit 53 of development apparatus (reception 
system) 50. Radio circuit 53 has an AGC (Automatic Gain 
Control) circuit and AFC (Automatic Frequency Control) 
circuit, and compensates for carrier frequency offset 
between transmission and reception and input level 

20 variations . The signal having undergone analog/digital 
conversion in analog BB processing section 52 is output 
to digital BB processing section 51. 

Digital BB processing section 51 is a portion to 
perform digital demodulation, despreading and other 

25 processingwhen de vel opment apparatus (reception system) 
50 is a CDMA (Code Division Multiple Access) reception 
apparatus, while being a portion to perform digital 
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demodulation, Fourier transform and other processing when 
the system 50 is an OFDM reception apparatus. The signal 
processed in digital BB processing 51 is input to error 
rate measurement instrument 70, and error rate 
5 measurement instrument 70 measures a channel error rate 
of the signal . 

Thus, in channel simulator 10, simulated multipath 
and fading variation on each path assumed to occur on 
transmission channels are added to radio signals obtained 

10 in transmission system 40 of the development apparatus, 
thus obtained signals are input to reception system 50 
of the development apparatus, and error rate 
characteristics of signals processed in reception system 
50 are measured, thereby evaluating transmission 

15 characteristics of transmission system 40 and reception 
system 50 . 

In recent years, as techniques enabling 
large-capacity data transmission, attention is being 
given to mul t i - an t enna techniques typified by MIMO (Multi 

20 Input Multi Output) and adaptive array antenna. For 
example, in a mu 1 1 i - an t e nna apparatus using the MIMO 
technique, a plurality of antennas is provided for the 
transmission and reception systems, each antenna of the 
transmission system transmits respective different data, 

25 and the reception system performs propagation path 
estimation and other processing on combined signals to 
separate, and restores a plurality of items of data. 
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In developing the multi-antenna apparatus, 
performance evaluations using a conventional channel 
simulator only enables inadequate evaluations . In other 
words, in the multi-antenna apparatus having M antennas 
5 on the transmission side and N antennas on the reception 
side, there exist MxN transmission channels, but the 
conventional channel simulator only allows measurements 
on a single channel. Even if the number of channels is 
increased, it is not extremely insufficient to evaluate 

10 the systems whose performance depend on spatial 
information such as transmission and reception antenna 
arrangements and radiation direction and the direction 
of arrival on each path. 

Further, in a channel simulator with the merely 

15 increased number of channels, in order to reproduce 
multi-channel using channel data collected in running 
experiments, it is necessary to collect data of all the 
channels and all the paths involved in the number of 
transmission and reception antennas and their 

20 arrangements of data collecting apparatuses used in a 
development apparatus and running experiments, and as 
well as the need of enormous memory to store data, running 
experiments should repeatedly be carried out whenever 
the number of transmission and reception antennas and/or 

25 their arrangement is changed. 

Disclosure of Invention 



It is an object of the present invention to provide 
a channel simulation method and channel simulator 
enabling MxN-channel transmission formed by a 
multi-antenna apparatus to be simulated with simplicity 
5 and excellence. 

The object is achieved by generating channel 
variations of all the channels based on information of 
arrangements of transmission and reception antennas . In 
the present invention, it is noted that a delay difference 

10 and phase difference of each path occur corresponding 
to antenna installation positions of transmission and 
reception antennas in each channel, and by varying only 
the delay difference and phase difference of each path 
in each channel, an MxN-channel variation model is 

15 simplified. 

Description of Drawings 

FIG.l is a block diagram illustrating a 
configuration of a conventional channel simulator; 
20 FIG. 2 is a view illustrating lxl-channel 

t r ansmi ssion; 

FIG. 3 is a view to explain paths; 

FIG. 4 (A) is a view illustrating a delay profile; 
FIG. 4(B) is a view illustrating an instantaneous 
25 variation; 

FIG. 4(C) is a view illustrating a short-term 
variation; 
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FIG. 4(D) is a view illustrating a long-terra 
variation ; 

FIG. 5 is a view to explain an elementary signal; 

FIG. 6 is a view illustrating a model with a single 
5 elementary signal; 

FIG. 7 is a view illustrating a case where a signal 
reflected diffusely by a sphere near a virtual antenna 
is received in 1 ine - o f - s igh t angle <f> as an elementary 
signal; 

10 FIG. 8 is a view illustrating propagation delays with 

a large line-of-sight angle; 

FIG. 9(A) is a view illustrating the direction of 
arrival of an elementary signal when the radius of a 
scattering sphere includes a reception antenna; 
15 FIG. 9(B) is another view illustrating the direction 

of arrival of an elementary signal when the radius of 
a scattering sphere includes a reception antenna; 

FIG. 10 is a view to explain a principle of generation 
of a standing wave under an environment with a large number 
20 of reflection waves; 

FIG. 11 is a view to explain a power density spectrum 
of variation in envelop amplitude due to Rayleigh fading; 

FIG. 12 is another view to explain a power density 
spectrum of variation in envelop amplitude due to Rayleigh 
25 fading; 

FIG. 13 is a view illustrating an MxN-channel 
transmission formed by a mu 1 t i - ant enna apparatus; 
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FIG. 14 (A) is a view to explain a difference in path 
distance caused by a distance between transmission and 
reception antennas, radiation angle, and the angle of 
arrival / 

5 FIG. 14(B) is another view to explain a difference 

in path distance caused by a distance between transmission 
and reception antennas, radiation angle, and the angle 
of arrival; 

FIG. 15 is a view illustrating an environment with 
10 signals arriving in all the directions; 

FIG. 16 is a view illustrating a model to add an 
instantaneous variation to each channel when multipath 
is not present; 

FIG. 17 is a view illustrating a model to add an 
15 instantaneous variation to each channel when multipath 
is present; 

FIG. 18 is a view illustrating a model to generate 
mutually correlated band-limited complex gaussian noises 
from MxNxP mutually independent band-limited complex 
20 gaussian noises; 

FIG. 19 is a block diagram illustrating a 
configuration to generate correlated instantaneous 
variations (two signals) proposed by Sasaoka; 

FIG. 20(A) is a view to explain a principle to form 
25 MxN-channel correlated instantaneous variations from 
lxl-channel instantaneous variations; 

FIG. 20(B) is another view to explain a principle 
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to form MxN-channel correlated instantaneous variations 
from lxl-channel instantaneous variations; 

FIG. 21 is a block diagram illustrating connections 
between a channel simulator and development apparatus 
5 according to an Embodiment of the present invention; 

FIG. 22 is a block diagram illustrating a 
configuration of a channel simulator of the Embodiment; 

FIG. 23 is a table showing descriptions of each 
parameter used in the Embodiment; 
10 FIG. 24 is a block diagram illustrating a 

configuration of a reference channel path control 
section; 

FIG. 25 is a block diagram illustrating a 
configuration of a channel processing section; 
15 FIG. 26 is a block diagram illustrating a 

configuration of a correlated gaussian noise generating 
section; 

FIG. 27 is a block diagram illustrating a 
configuration of a reference channel path control 
20 section; 

FIG. 28 is a block diagram illustrating a 
configuration of a channel processing section; 

FIG. 29 is a block diagram illustrating a 
configuration of a correlated gaussian noise generating 
25 section; 

FIG. 30 is a block diagram illustrating a 
configuration of a fading adding section; 
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FIG. 31 is a block diagram illustrating a 
configuration of a transmission analog adjusting section; 

FIG.32 is ablock diagram illustrating a pseudo power 
amplifier (PA); and 
5 FIG. 33 is a block diagram illustrating a 

configuration of a reception analog adjusting section. 

Best Mode for Carrying Out the Invention 

In order to accurately simulate transmission 

10 channels in a multi-antenna apparatus having M 
transmission antennas and N reception antennas, different 
channel variations should be added to MxN transmission 
channels. However, when channel variations of MxN 
channels are simulated by simply adding channel variation 

15 parameters for each channel, enormous amounts of 
parameters and computations are required and an apparatus 
configuration becomes complicated. 

The inventor of the present invention thought that 
if channel models are simplified in forming the channel 

20 models of MxN channels in a multi-antenna apparatus, the 
number of parameters and computation amounts would be 
reduced, and as a result , an apparatus configuration would 
also be simplified relatively, and has reached the present 
invention . 

25 It is a gist of the present invention generating 

channel variations on all the channels based on 
information of arrangements of transmission and reception 
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antennas. In the present invention, it is noted that 
a delay difference and phase difference of each path are 
generated corresponding to antenna installation 
positions of transmission and reception antennas in each 
5 channel, and channel variation models of MxN channels 
are simplified by varying the delay difference and phase 
difference of each path in the channels. Further, in 
the present invention, on the assumption that a correlated 
instantaneous variation is multiplexed on a signal of 
10 each path, a channel variation model is generated such 
that the correlated instantaneous variation is added for 
each path . 

By this means, it is possible to form channel 
variation models corresponding to a plurality of channels 

15 with a relatively small amount of calculations using 
approximate multipath channels obtained from existing 
single channel measurement data calculated in RayTrace 
simulation, actual running experiments, and so on. 

Further, in the present invention, as a method of 

20 generating instantaneous variations (correlated complex 
gaussian noises) correlated between channels or paths, 
five following methods are proposed: 
® Eigenvalue transformation (spatio-temporal) 

(5) Eigenvalue transformation (spatio) 

25 (3) Cholesky factorization (spatio-temporal) 

® Cholesky factorization (spatio) 
(5) Expanded Sasaoka method (spatio-temporal) 
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Among these four methods, items 1 to 4 are to devise 
a method of obtaining a transformation matrix A in 
calculating a correlated instantaneous variation to add 
to each channel from mutually independent band-limited 
5 complex gaussian noises using the transformation matrix 
A. Item 5 is to devise to enable generation of correlated 
instantaneous variations on MxN channels by expanding 
the Sasaoka method of generating correlated instantaneous 
variations proposed on 1x2 channels. 

10 (1) Principles in this Embodiment 

Prior to descriptions on a constitution of the 
Embodiment, the principles in the Embodiment will be 
described first. The inventor of the present invention 
first considered different points and similar points 

15 between Ixl-channel transmission and MxN-channel 
transmission. Further, in order to expand a lxl channel 
transmission model to an MxN-channel transmission model 
with simplicity as possible, detailed studied were 
conducted on how to expand a short-term variation and 

20 instantaneous variation. The studies will sequentially 
be described below. 

(1-1) lxl-channel transmission 

FIG. 2 illustrates lx 1 - channe 1 t ran smi s s i on ea ch wa y 
between a transmission antenna and a reception antenna. 
25 Hereinafter, a channel between one-to-one transmission 
and reception antennas is referred to as 1 channel . FIG. 3 
illustrates paths . While the channel is shown by 
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straight line in FIG. 2, signals are actually reflected 
and diffracted in space, thereby passed through various 
paths ((D-@ in FIG . 3) , and received in a receiver. Then, 
propagation delays vary with path distances, and a delay 
5 profile is drawn with propagation delay time on the 
horizontal axis and reception power on the vertical axis, 
as shown in FIG. 4 (A). Signals arriving with different 
delays must have been passed through different 
propagation paths, and the propagation paths are simply 

10 referred to as paths. 

Each path is defined by transmittal coefficient 
(complex) indicative of extents of delay, gain (actually, 
attenuation) and phase shift to which a signal passed 
through the path undergoes. Measuring the delay profile 

15 enables acquisition of the approximate number of paths 
constituting the channel, and of scales of the delay and 
gain that each of the paths provides. A phase on each 
path varies with traveling speed and the angle of arrival 
of a signal with respect to the traveling direction. 

20 FIGs.4(B) to 4(D) illustrate gain variations of a 

path (note that the horizontal axis represents the 
distance, instead of time). Gain variations are 
classified into a long-term variation (distance 
variation) dependent on the distance from the 

25 transmission antenna and di r e ct i ona 1 i t i e s of 
transmission and reception antennas, short-term 
variation (shadowing) due to e f fects of shields by obj ects 
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on the ground, and instantaneous variation due to 
mult i- wave multiplexing . 

The propagation distance and propagation delay are 
in a proportional relationship, and therefore, the 
5 long-term variation has almost the same shape as that 
of the delay profile. When a terminal travels in mobile 
communications, the propagation distance (or propagation 
delay) of each path varies, and the reception level also 
varies, but the rate of the variation is the slowest (very 

10 slow) among the other variations. The long-term 
variation was mode led with the Okumura curve (Hata method ) 
which was made by analyzing statistically a lot of running 
experiment data and has been used widely. Further, in 
recent years, the Sakagami method has been used which 

15 is modified by adding a use frequency band and parameters 
of ground objects. 

The short-term variation is a gain variation caused 
by each path being shield or appearing by/from a building 
or the like (a path may be shield by a person walking 

20 around in wireless LAN) . With respect to the variation 
rate, there is no theoretical equation, but in general, 
the rate is assumed to be 1Hz or less. Actually, the 
short-term variation should be determined from a cause 
of occurrence in relation to objects on the ground and 

25 traveling speed. For example, when a terminal travels 
at speed of 30Km/h in a building street with a building 
width of 30m, it is considered that a variation occurs 
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in a cycle of 3.6 seconds, and surely, such a frequency 
is estimated to be 1Hz or less in many cases. The gain 
variation due to the short-term variation follows a 
logarithm normal distribution, and is modeled such that 
5 the gain varies concurrently in a band (in the above case, 
the band ranges from 0 to 1/3.6 [Hz]). 

The instantaneous variation is a variation 
occurring when several elementary signals are multiplexed . 
On a path looking like a single path on the delay profile, 

10 a plurality of signals (such that the amplitude and phases 
are not in complete agreement) is passed actually. The 
multiplexed signal thus seeming to pass on a single path 
on the delay profile is referred to as an elementary signal, 
and the amplitude and phase varies (which is regarded 

15 as a single signal-wave being passed through a path 
providing variations in gain and phase) . The 
instantaneous variation is explained by the Doppler 
effect, and varies at speed of about several hertz to 
1kHz as described later. 

20 Properties of an elementary signal causing the 

instantaneous variation are characterized by the angle 
of arrival 9 and line-of -sight angle (j> (variation width 
of the angle of arrival 9 ) . In FIG. 5, attention is 
directed toward elementary signal (3) in FIG. 4 (A). Paths 

25 are widely divided in space in FIG. 4 (A), and therefore, 
the receiver receives signals undergoing completely 
different variations. Accordingly, for example, 
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elementary signal ® can be regarded as a signal «in 5 

from a virtual transmission antenna on an extended 1». 

rt f arrival (on the assumption that the 
in the direction of arrival iu 

•™ at soeed v upwardly as viewed in the 
receiver is moving at speea 

3 figure) . 

e~ of a sinqle elementary signal. 
FIG. 6 shows a case of a smgx« 

In this case, the elements signal is not a multiplexed 
signal, and is received without any variations rn 
amplitude and phase except coppler shift doe to traveling . 
0 such a case hardiy occurs in mobile communications, but 
is so.eii.es used as a channel »oei. In addition, an 
amount of Doppler shift of elementary signal © rs 
represented by f.cos 8 "here f.-vxfc/c 

15 by . sphere near a virtual antenna is received in a 
U».-of-.i,ht angle * (variation width of the angle of 
atrival 6, as an elementary signal. In this case, the 
elementary signal varies in amplitude and phase, but the 
ang le of arrival B can be measured stably because the 

• it a nH the delay difference 

20 li„e-of-sight angle * rs small, and the y 

i. thus small also. However, as the 1 ine-of -s ight angle 
* is increased (the radius of the scattering sphere near 
th . virtual antenna is increased, , the elementary signal 
includes signals with considerably large delay 
25 di££ ere„cesasshow„rnr IG .8. as , re suit , the amplitude 
and phase both vary largely, and it becomes difficult 
to measure the angle of arrival 6 . 
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When the radius of the scattering sphere is further 
increased and starts including the reception antenna as 
shown in FIG. 9 (A), the elementary signal seems to come 
from in every direction, and in this case, should be 
5 divided into elementary signals of a plurality of paths 
on the delay profile. However, even if the elementary 
signal is divided on the delay profile, divided signals 
also contain a number of signals which have same 
propagation delays but are passed through spatially 

10 completely different paths, and therefore, each of 
divided elementary signals seems to come from every 
direction. In other words, an image as shown in FIG. 9(B) 
is obtained (like the figure with the inverse transmission 
direction to FIG. 7). As a matter of course, it is not 

15 possible to measure the angle of direction of each 
elementary signal (such measurement is meaningless) . 

In an environment with a large number of reflection 
waves as shown in FIG. 10, a standing wave occurs, and 
with respect to the reception level, spatially strong 

20 points and weak points appear repeatedly. This is the 
reason of spatial distribution of the instantaneous 
variation. In particular, the instantaneous variation 
occurring when independent signals of almost the same 
level arrive in all the directions is referred to as 

25 Rayleighfading, anditis known that the envelop amplitude 
variation conforms to the Rayleigh distribution and that 
the phase variation conforms to the uniform distribution. 
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FIGs.ll and 12 illustrate power density spectra of 
the envelop amplitude variation due to Rayleigh fading. 
In the case where signals arrive in directions CD,®, 
(3) and ® with respect to the traveling direction v as 
5 shown in FIG. 11, the signal coming in direction ® seems 
to have the highest frequency, and reversely, the signal 
coming in direction ® seems to have the lowest frequency. 
This maximum frequency deviation is referred to as maximum 
Doppler frequency f D . The maximum Doppler frequency f D 

10 can be calculated as the number of standing waves contained 
in a traveling distance in a second where the standing 
waves are repeated in a cycle of wavelength, and generally, 
ranges from about several hertz to 1kHz (when the carrier 
frequency is 2GHz and the traveling speed is lOOkrn/h, 

15 f D is a multiplication of these numerals, 200Hz). 
Similarly, a frequency deviation of f D • cos 6 is imposed 
on the signal © coming at angle 8 toward the traveling 
direction v, while the signal (3) with 6 of 90 degrees 
does not undergo the deviation. 

20 In addition, for example, in the instantaneous 

variation such that elementary signal ® arrives only in 
shaded directions (the angle of arrival 0= 180°, the 
line-of-sight angle (/>= 80°) as shown in FIG. 11, when five 
or more signals are contained, the variation is regarded 

25 as Rayleigh fading. The power density spectrum is only 
a shaded potion as shown in FIG. 12. 

According to the foregoing, it is possible to define 
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a long-term delay profile, short-term delay profile and 
instantaneous delay profile. Each path on the long-term 
delay profile undergoes the long-term variation. Its 
delay and gain are determined according to the condition 
5 of objects on the ground, traveling speed and direction 
and the angle of arrival, and vary slowly. Each path 
on the short-term delay profile undergoes the short-term 
variation (shadowing) in addition to the long-term 
variation. By this means, the gain of each path varies 

10 at a rate of 1Hz or less in the logarithm normal 
distribution independent of each path. 

Each path on the instantaneous delay profile 
undergoes the instantaneous variation in addition to the 
long-term variation and short-term variation. By this 

15 means, the gain and phase of each path are subjected to 
Rayleigh fading (with the gain of the Rayleigh 
distribution and the phase of the uniform distribution) 
independent of each path. The variation rate is 
determined by the carrier frequency, traveling speed, 

20 the angle of arrival and 1 ine - o f - s i gh t angle, and ranges 
from several to several hundred hertz. 

In addition, the amplitude on the delay profile 
represents reception power of an elementary signal coming 
from each path, and does not have the gain nor phase 

25 (accurately, which is complex amplitude of complex 
impulse response on each channel) . Conversely, 
"'representation of power of each path" is not appropriate , 
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but may be used in the scope without causing 
misunderstanding according to convention. 

(1-2) Expansion to MxN-channel transmission 
(1-2-1) Expansion of the short-term variation 
5 FIG. 13 illustrates MxN-channel transmission formed 

by a mul t i -an tenna apparatus having M transmission 
antennas and N reception antennas. 

The inventor of the present invention thought that 
MxN channels are similar to one another. In other words, 
10 unless the arrangements of transmission and reception 
antennas are spread out in an area of several meters square, 
assuming that an interval of the short-term variation 
is about ten and several meters, it should be regarded 
that not only long-term delay profiles but also short-term 
15 delay profiles respectively observed actually in the 
reception antennas are almost equal to one another. 

Then, considerations reached that differences 
between channels on the short-term delay profile are only 
in propagation delay and carrier phase of each path due 
20 to a difference in path distance resulting from the 
arrangements of transmission and reception antennas. 

FIG. 14 (A) illustrates comparison between paths from 
two transmission antennas to a single reception antenna. 
When distance d T between elements is small enough, paths 
25 from vicinities of the antennas to the reception antenna 
can be regarded as being common, a difference in path 
distance is thus d T ■ cos0 T , and corresponding to the 
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difference, the channel has differences in path delay 
and phase (that is a carrier phase, but can be referred 
to as a path phase) . 

Reception antennas as shown in FIG. 14(B) have 
5 similar phenomenon. However, note that the radiation 
angle 6 T and the angle of arrival 9 R are defined for 
each path. In the case where radiation angles and angles 
of arrival range in all the directions as shown in FIG. 9, 
differences in path distance are reversed with angles, 

10 and on average, can be considered as being not present. 
It is the same as in a path with unknown radiation angle 
and angle of arrival. 

Based on the aforementioned considerations, the 
inventor of the present invention has reached a conclusion 

15 that if the radiation angle and angle of arrival are 
obtained in each path on the short-term delay profile 
of either one of the channels, short-term delay profiles 
of the other channels can be calculated from arrangements 
of transmission and reception antenna elements. 

20 Accordingly, in the present invention, using 

existing one-channel transmission measurement data 
(including short-term delay profiles) obtained in 
Ray-Trace simulation, actual running experiments and so 
on, short-term delay prof iles (variations in delay, gain 

25 and phase of all the paths) of MxN channels are formed 
by calculations from the one-channel transmission 
measurement data and arrangements of transmission and 
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reception antenna elements, thereby forming channel 
variation models of MxN channels . It is thereby possible 
to generate information of paths on all the channels with 
simplicity and accuracy from information of each path 
5 on the reference channel. 

(1-2-2) Expansion of instantaneous variation 
Considered next is an instantaneous variation on 
each channel. The inventor of the present invention 
studied on how instantaneous variations of channels and 

10 paths are similar to one another. 

In the case of the instantaneous variation, for 
example, with respect to standing waves with a band of 
5GHz, on average, a distance of 3cm ( ha 1 f- wave length ) 
contains a single wave . Therefore, unlike the short-term 

15 variation, instantaneous variations are not equal on all 
the antennas . However, it may be said that at the moment 
when some antenna undergoes a position variation, an 
adjacent antenna generally undergoes a position variation 
also. The similar phenomenon may occur in time. The 

20 instantaneous variation varies at a rate of 1kHz or less, 
and when a positive variation is imposed at a moment, 
it is likely that the positive variation lasts 0 . 1ms later. 
The former is quantitatively represented by spatio 
correlation function, while the latter is quantitatively 

25 represented by temporal correlation function, and with 
respect to Rayleigh fading where signals come from all 
the directions, these functions are led to following 
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equations : 
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Spatial correlation p[d)=x .* (t)-xj{t) «J 0 (2*<*M) 

function: (-j) 

Temporal correlation r) = — (t) . x . [t+T) = j 0 (2nf D r) 
function: 



Temporal correlation function: 

,-n x(t) and Xj (t) respectively 
In equation ID f X ^ L) 

• H .ianals (complex baseband signals) 
represent received signals i 

.q d is a distance between the 
of ith and jth antennas, d is a 

f = nath X is wavelength, 

thet . 

and f. is the ooppier fluency. Further, 

represents a conjugate »,le. n-e, ana 1 0 represents 

Bessel function. 

th e case of conventional equations on 
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pet£ or m ance o £ antenna diversity, h y settin, a distance 
b et»een antennas so that the antennas are 
spaced, tor e*a m ple, hy h»lf-wa,elen g th such that the 
apatio correction value is stall, evaluations have heen 
carr ied out on the assumption that received si.nals are 
yarded as hein, not correlated with one another. rhrs 
ls ba sed on the spatio correlation .unction o £ e q uatro„ 

that the time correlation is 
(1), but it is ignored that 

actually present. 

in order to obtain more accurate 
Accordingly, m oraer 
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results, there is needed a theory enabling concurrent 
analysis of temporal correlation and spatio correlation. 
With respect to this issue, Sasaoka derived that the 
spatio-temporal correlation function becomes the 
5 following equation under circumstances as shown in FIG. 15 
(Sasaoka, "New Generation Method of Mutually Correlated 
Multipath Fading Waves " , IEICE, Transaction, 88/6, 
Vol, J71-B NO. 6) . Herein, in FIG. 15, </> represents an 
arrangement angle of an antenna with respect to the 
10 traveling direction. 

S patio-temporal 
correlation function: 

P(d,T) = J 0 [V ( 27rd r /tf+(27i;f D T r ) 2 ] 

(2) 

where d r = d-siny/ , r r = r-(d / Af D )-cosy/ 



In addition, FIG. 15 shows an environment where 
15 signals come from all the directions. In the case where 
signals come from only some direction, an instantaneous 
variation may be used which has only part of U-shaped 
powe r density spectrum (FIG. 12) o f Ra yl e i gh fading . This 
is because, while it is obvious to a sinusoidal wave, 
20 correlation between waves with different frequencies is 
"0" spatially and temporally, and therefore, equation 
(2) holds independently, irrespective of frequency 
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components (i.e. the direction of arrival) . 

Considered next is the correlation between 
multipath (accurately, multipath that is identified by 
delay time) . For example, as shown in FIG. 5, in the case 
5 where an elementary signal passed through each path is 
a single signal, since the elementary signal is a complete 
beam and the instantaneous variation does not occur, the 
short-term delay profile becomes the instantaneous delay 
profile. In other words, the correlation does not occur 
10 between the multipath by beam, which corresponds to the 
case where the line-of -sight angle is w 0", and therefore 
the correlation is "0" because instantaneous variation 
spectra are not in agreement unless the 1 ine-o f - s ight 
angles are in agreement. 
15 In contrast thereto, in the case where the 

line-of -s ight angle is extremely large as shown in FIG. 9, 
i.e. an elementary signal includes signals arriving from 
all directions, it is considered that in paths with large 
delay differences, spatial path differences between the 
20 paths become large and the correlation between the paths 
should be almost zero, while in paths with small delay 
differences, spatial path differences are small and 
similar instantaneous variations are imposed on the paths 
(i.e. there is the correlation) . In such a case, it is 
25 considered that equation (2) can also be applied to paths. 

In the case of a path with the instantaneous 
variation of an almost unified elementary signal such 
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that the line-of - s ight angle is in a certain range as 
shown in FIG. 7 , it is considered obtaining an intermediate 
result between the two paths described above. In other 
words, following descriptions are considered. 
5 (a) In the case of a few overlaps of the angle of 

arrival and line-of -sight angle in paths, the paths are 
obviously different in spatial route, and have small 
correlation (where the extreme case indicates "beam" in 
FIG . 5 ) . 

10 (b) In the case where overlaps of the angle of arrival 

and line-of -sight angle are large and the difference in 
delay time is also large in paths, the paths are also 
different in spatial route, and have small correlation 
(which corresponds to the case of FIG. 9) . 

15 (c) The correlation between paths received in 

different antennas largely depends on the distance 
between the antennas. 

As a result, with respect to the correlation between 
paths, (a) can be implemented by performing band 

20 limitation on instantaneous variations corresponding to 
the angle of arrival and the 1 ine - o f - s i gh t angle, because 
different paths make the correlation small. (b) and (c) 
can be implemented by applying the spatio-temporal 
correlation function of eq. (2) . 

25 In the view of the foregoing, the MxN-channel 

transmission as shown in FIG. 13 without multipath is 
represented by MN channels as shown in FIG. 16. Herein, 
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the gain of the short-term variation of each channel is 
assumed to be equal, and thus, FIG. 16 shows only the 
instantaneous variation of each channel. Further, in 
order to eliminate effects of transmission data in 
5 observing the similarity degree between channels, 
transmission data is fixed to "1". Then, the variation 
is added by multiplying a signal of each channel by an 
instantaneous variation (complex gaussian noise) with 
spatio correlation based on a distance between antennas. 

10 In a case of p-multipath, simulation is expanded 

as shown in FIG. 17, and a signal of each path is multiplied 
by an instantaneous variation with spatio-temporal 
correlation of eq. (2) based on the distance between 
antennas and delay difference. In other words, in a 

15 configuration as shown in FIG. 17, instantaneous 
variations are added using band-limited complex gaussian 
noises in consideration of the angle of arrival and the 
line-of -sight angle so that each path complies with the 
spatio-temporal correlation function of (2) . 

20 (1-3) Methodof generating correlated instantaneous 

variations 

The issue is how to generate desired correlated 
complex gaussian noises. As described above, as a method 
of generating correlated instantaneous variations, 
25 proposed in this Embodiment are correlation matrix 
methods (eigenvalue transformation method and Cholesky 
factorization method) and expanded Sasaoka method. 
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FIG. 18 illustrates a method of generating MxNxP (P: 
the number of paths) mutually correlated band-limited 
complex gaussian noises ( co r r e la t ed gaus s ian no i s e s ) from 
MxNxP mutual ly independent band-limited complex gaussian 
noises, using transformation matrix A (with MNP rows and 
MNP columns). The issue is to obtain transformation 
matrix A that provides desired path correlation. 

Assuming that outputs of paths are Y=(yi, 
y2/ • • • / yMNp) T / the path-correlation matrix is as indicated 
in the following equation. In addition, for simplicity, 
numeral subscripts are indicated by serial numbers below. 
Further, in the following equation, * in upper subscript 
represents a conjugate complex number, H represents 
conjugate complex transpose, and E() represents an 
ensemble mean. 



Each of (MNP) 2 elements represents correlation 
between paths, and is obtained by calculating a 
spatio-temporal correlation value of eq. (2) from a 
difference in path distance and a difference in 
propagation delay obtained from transmission and 




(3) 
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reception antenna arrangements and the radiation angle 
and the angle of arrival of a signal, and thus, desired 
path correlation matrix S Y y is obtained. 

In other words, it is required to generate Y such 
5 that 2 yy =E(YY H ), and to generate such Y, the issue is 
how to determine transformation matrix A in FIG. 18. This 
Embodiment proposes the method based on eigenvalue 
transformation, and the method based on Cholesky 
factorization . 

10 (1-3-1) Eigenvalue transformation method 

In FIG. 18, such a relationship that Y = AX holds 
between complex gaussian noises. Herein, if an inverse 
matrix of transformation matrix A exists, the following 
equation holds: 

15 X = A _1 Y ... ( 4 ) 

Eq. (4) represents obtaining (MNPxl) signal vector 
X without mutual correlation by multiplying (MNPxl) 
signal vector Y with correlation as in eq. (3) by matrix 
A" 1 , and this relationship is known as eigenvalue 

20 transformation (or KL transformation) . In eigenvalue 
transformation, used as A -1 is MNP (MNPxl) unit 

e igenvect o r s ei , e 2 , , e MNP , in a row, of desired path 

correlation matrix 2 Y y, and transformation matrix A is 
indicated in the following equation: 
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A =(<?! e 2 




a-(a-')" .(a-')"- 



(5) 



By thus obtaining transformation matrix A by the 
eigenvalue transformation method, it is possible to use 
a matrix with a small number of elements in obtaining 
correlated instantaneous variations from mutually 
independent instantaneous variations, and thus, the 
correlated instantaneous variations can be obtained with 
a small amount of calculations. 

In addition, described herein is the case of forming 
mutually correlated instantaneous variations that are 
correlated between channels and between paths by forming 
MxNx (the number of paths) correlated instantaneous 
variations correlated between paths by performing matrix 
operationprocessingonMxNx (the number of paths ) mutual ly 
independent instantaneous variations using 

transformation matrix A . However, thepresent invention 
is applicable without being limited to such a case. 

For example, correlated instantaneous variations 
corresponding to MxN channels that are correlated between 
all the channels may be formed by forming correlated 
instantaneous variations corresponding to MxN channels 
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correlated with channels by performing the matrix 
operation processing on a plurality of mutually 
independent instantaneous variations corresponding to 
MxN channels using transformation matrix A. The 
5 foregoing is the same as in following descriptions. 
(1-3-2) Cholesky factorization method 
Desired path- correlation matrix 2 yy can toe 
subjected to Cholesky factorization as in the following 
equation . 

10 E yy =L H L - - - ( 6) 

where L is an (MNP x MNP) lower triangular matrix. 

Transformation matrix A in FIG. 18 is obtained using 
obtained lower triangular matrix L in the following 
15 equation . 

A = L H ... (7) 

This is because path correlation matrix E(Y H Y) is 
indicated as in the following equation, since Y=AX=L H X 
from FIG. 18. 

2 0 E ( YY H ) =E [ (L H X) ( L H X ) H ] =E ( L H XX H L) =L H E (LXX H ) L = L H L= 2 Y y 

... (8) 

Herein, since each element of vector X is an independent 
gaussian variable, correlation matrix E(XX H )=I is used. 

With respect to Cholesky factorization, an 
25 approximate algorithm enabling large reductions in 
computation amount has been disclosed recently 
(H.R.Karimi etc. : "A Novel and Efficient Solution to 
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Block-Based Joint - Detect ion Using Approximate Cholesky 
Factorization", P IMRC ' 9 8 . p , 1 3 4 0 - 1 3 4 5 , 1998), and it is 
effective using such an algorithm. 

By thus obtaining transformation matrix A using 
5 Cholesky factorization, it is possible to use a lower 
triangular matrix obtained by Cholesky factorization in 
obtaining correlated instantaneous variations from 
mutually independent instantaneous variations, and it 
is thereby possible to obtain the correlated 
10 instantaneous variations with a small amount of 
calculations . 

(1-3-3) Expanded Sasaoka method 

Proposed next is the method (hereinafter, referred 
to as an expanded Sasaoka method) of expanding the Sasaoka 

15 method to MxN-channel transmission, aside from obtaining 
correlated instantaneous variations from matrix A. 

The method proposed by Sasaoka is first explained 
briefly. FIG. 19 illustrates a block diagram to generate 
correlated instantaneous variations (two signals) 

20 proposed by Sasaoka . Inaddition, s ys terns and pa ramete rs 
assumed in FIG. 19 are the same as those as shown in FIG. 15. 
With respect to FIG. 19, Doppler filters 102 and 105 perform 
spectral shaping on two-system white gaussian noises 
generated in white gaussian noise generating sections 

25 101 and 104 to power density spectra of Rayleigh fading 
when an elementary signal includes signals arriving in 
all the directions, respectively (when the angle of 
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arrival and the 1 ine - o f - s i gh t angle of an elementary 
signal are know, corresponding to which, the bandwidth 
i s narrowed ) . 

In order to cause the noises to be correlated with 
5 each other, the noises are passed through two types of 
filters, 103 and 106, respectively. In the conventional 
generation method, filters 103 and 106 have spatio 
correlation value p and V~ (1- p 2 ) , and could not express 
temporal correlation. In order to also provide temporal 

10 correlation, Sasaoka replaced the filters. Filter 
characteristics H(f ) and G(f) are determined by distance 
d between two reception antenna elements, angle (/> of 
antenna arrangement to the traveling direction, carrier 
wavelength X , and maximum Doppler frequency f D , and since 

15 gains are in sin-cos relation, the power density spectrum 
of an instantaneous variation of elementary signal 2 to 
input to complex multiplier 111 has the same shape (shape 
of Doppler filter) as that of elementary signal 1. 
Finally, corresponding to a sign of cos 0 , either delay 

2,0 circuit 108 or 109 provides one of the elementary signals 
with delay. In addition, an output from delay circuit 
108 is valid when cos 0 2^0, while an output from delay 
circuit 109 is valid when cos <J> < 0. 

Filters 103 and 106 in FIG. 19 have complicated 

25 characteristics (that vary with conditions) , and Sasaoka 
proposed to implement the filters by performing weighted 
addition of multitone with different frequencies. 
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Then, by spectral shaping as shown in FIG. 19, filter 
characteristics H(f) and G(f) of filters 103 and 106 are 
obtained to generate correlated instantaneous variations 
meeting the spatio-temporal correlation function of 
eq. (2) (with respect to correlated instantaneous 
variations prior to Sasaoka, these filers having spatio 
correlation coefficients have been proposed, but it is 
descried that excellent approximation is not obtained 

unless d/ /LCI). 

In FIG. 19, two output signals of the Doppler filters 
are assumed to be xi (t) and x 2 (t) . In order not to depend 
on transmission signals because of study on correlation 
between channels, it is further assumed that inputs of 
two elementary signals are " 1 " , and that outputs of the 
elementary signals provided with instantaneous 
variations are yi(t) and y 2 (t) (i.e. which are the 
correlated instantaneous variations) . Each power 
density spectrum (expression in frequency of the 
auto-correlation function) and cross power density 
spectrum (expression in frequency of cross-correlation 
function) are indicated by following equations: 

Sylyl ( f ) = Sxlxl ( f ) 

= S(f) : Doppler Spectrum 

s y2y2 (f)= |h (f ) l 2 s xlxl (f ) + |g (f ) | 2 s x2 x2 (f ) 
= { |h (f ) l 2 + |g (f ) l 2 }s (f ) 

Syly2 (f ) = H* (f ) S X 1 X 1 (f) 

= H* (f ) Sxixi (f) ... (9) 
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The correlated instantaneous variations should have 
power density spectral with the same shape, and therefore, 
meet the following equation: 
5 |H(f ) | 2 +|G(f ) 1 2 =1 ... (10) 
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Further, cross power density spectrum S y i y 2(f) is 
obtained by subjecting eq. (2) to Fourier Transform, and 
the following equation holds: 



Hif) =cos\ 



27id'sin 




(11) 



Then, with respect to G(f), the following equation 
is obtained from equations (10) and (11) : 



G if) =sin 1 



27id'sin 



(12) 

15 

In eq. (12) , a positive or negative sign appears in 
obtaining a square root, and since either sign does not 
affect a correlation value of instantaneous variations, 
the positive sign is selected. 
20 Considered below is a method of applying the method 

of generating correlated instantaneous variations as 
shown in FIG. 19 to MxN-channel transmission. FIGs . 20 (A) 
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and 20(B) are illustrations of the principle of forming 
MxN-channel correlated instantaneous variations from 
lxl-channel instantaneous variation. 

FIG. 20 (A) illustrates lxN channels and lxM channels, 
5 and FIG. 20(B) illustrates conversion from lxM channels 
into Mxl channels. Considered first is the case of 
receiving signals from a single transmission antenna at 
N reception antennas as shown in FIG. 20(A) (in this case, 
received signals likely have the correlation) . It is 

10 herein assumed that complex impulse response is known 
on a transmission channel between the transmission 
antenna and reception antenna 1 (hereinafter, such a 
channel is referred to as a 1-1 channel), and that the 
short-term delay profile of each channel has been 

15 calculated. 

With respect to instantaneous variations, the 
method in FIG. 19 is applied to each path of the 1-1 channel 
and each corresponding path of a 1-2 channel . Thereafter, 
in the same way, correlated instantaneous variations of 

20 a 1-3 channel, , 1-N channel are generated while 

varying the distance between antennas. 

Shown on the left of FIG.20B is lxM-channel 
transmission from the right to the left as viewed in the 
figure. Assuming that the transmission and reception 

25 antennas are in the same positions, the 1-1 channel has 
the same complex impulse response as that of the 1-1 
channel of FIG. 20(A) from the reversibility of channel. 
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Accordingly, it is possible to generate all correlated 
instantaneous variations on the left of FIG. 20(B) as in 
FIG. 20 (A) . Then, by changing the signal direction using 
the reversibility of channel again, correlated 
5 instantaneous variations of all the channels on the right 
of FIG. 20(B) are obtained (in other words, even if a 
transmission antenna is different, the channel has 
correlation when the reception antenna is just one) . 

By applying the left-hand part of FIG. 20(B) to 

10 reception antennas of FIG. 20 (A), it is possible to 
generate all the correlated instantaneous variations on 
MxN channels in transmission. However, it is necessary 
to note following respects. 

1 . Use independent white gaussian noises in paths 

15 that are divided and identified on the delay profile. 
This is because different paths are considered as 
spatially distinguished paths. 

2. While the 1-1 channel is used as a reference in 
the foregoing, any channel is available as a reference 

20 channel. This is because changing the reference channel 
changes propagation delays and phases, but does not 
provides changes in relative value. In addition, when 
antenna elements are circularly arranged, it is possible 
to set a reference channel in a center position where 

25 an antenna is not present actually. 

3. The instantaneous variation of each channel is 
defined by correlation with the reference channel and 
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generated, and for example, in FIG. 20 (A) , the correlation 
between the 1-2 channel and 1-3 channel is not ensured. 
In other words, in consideration of the fact that the 
correlation is a cosine function value between two data 
5 vectors, this expansion method is not appropriate. 
( 1- 4 ) Summary 

As proposed in the foregoing, it is possible to 
calculate delays, phases and correlated instantaneous 
variations of paths on all the channels from the 
10 installation information of transmission and reception 
antennas, and generate all the channel variation models 
on MxN channels. It is thereby possible to use existing 
single-channel transmission measurement data obtained 
in Ray-Trace simulation, actual running experiments or 
15 the like as MxN-channel transmission measurement data. 
(2) A configuration of this Embodiment 
(2-1) Entire configuration 

FIG. 2 1 illustrates connections of channel simulator 
120 according to this Embodiment and development 
20 apparatuses 40 and 50, with the same sections as in FIG.l 
assigned the same reference numerals as in FIG.l. In 
addition, descriptions are omitted herein on portions 
already explained using FIG.l. 

Channel simulator 120 simulates channels of 
25 development apparatuses 40 and 50 having a multi-antenna 
configuration, and thereby enables evaluations of channel 
characteristics of development apparatuses 40 and 50. 
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Channel simulator 120 is capable of receiving as 
its inputs digital baseband signal DB from digital BB 
processing section 41 of transmission system 40, analog 
baseband signal AB from analog BB processing section 42, 
5 and radio signal RF from radio circuit 43. An output 
of simulator 120 is selectively output to digital BB 
processing section 51, analog BB processing section 52 
or radio circuit 53 of reception system 50 corresponding 
to the operation of switches SW3 and SW4 . 

10 By this means, it is possible to directly input 

digital baseband signal DB from digital BB processing 
section 41 to channel simulator 120 even if development 
of radio circuit 43 and analog BB processing section 42 
has not been completed up to an operable state, and it 

15 is thereby possible to evaluate channel characteristics 
of digital BB processing sections 41 and 51 independently. 

As a result, without waiting for completion of radio 
circuits 43 and 53 (particularly, radio circuit 53 of 
reception system 50), it is possible to check the operation 

20 of digital baseband processing sections 41 and 51 that 
perform principal processing, whereby efficiency in 
development can be improved. 

FIG. 22 illustrates a configuration of channel 
simulator 120. In channel simulator 120, interface 

25 section 122 receives radio signal RFin from radio circuit 
43, analog baseband signal ABin from analog BB processing 
section 42 or digital baseband signal DBin from digital 
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BB processing section 41. More specifically, radio 
signals RFin or analog baseband signals ABin 
corresponding to the number (M) of transmission antennas 
are input to analog circuit 123, converted into digital 
5 baseband signals in analog circuit 123, and output. 
Switch SW10 selects either input digital baseband signal 
DBin or the digital baseband signal converted in analog 
circuit 123 to output to transmission analog adjusting 
section 124. In addition, the baseband signal is 

10 comprised of an I signal and Q signal, thereby forming 
2xM signals, and thus, indicated by 2M in the figure. 

In other words, in circuits subsequent to 
transmission analog adjusting section 124, subjects of 
processing are M digital baseband signals. More 

15 specifically, M transmission analog adjusting sections 
124 are provided corresponding to the number of digital 
baseband signals, and compensate M digital baseband 
signals for changes in transmission characteristics 
occurring due to fluctuations in performance of M analog 

20 BB processing sections 42, radio circuits 43 and analog 
circuits 123 of development apparatus (transmission 
system) 40. Specific configurations of transmission 
analog adjusting sections 124 will be described later. 

Switch 12 5 as signal duplication means makes N copies 

25 of each of M digital baseband signals, thereby forms MxN 
digital baseband signals, and outputs the signals to MxN 
channel processing sections, 126-1 to 126-MN. Each of 
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channel processing sections 126-1 to 126-MN receives 
information such as channel model information of a 
reference channel and transmission and reception antenna 
arrangement information formed in reference channel path 
5 control section 127, constructs a channel model of an 
assigned channel, and then, provides the digital baseband 
signal of the assigned channel with short-term complex 
impulse response and correlated instantaneous variation 
for the assigned channel corresponding to the constructed 

10 channel model in a complex multiplier. Specific 
configurations of channel processing sections 126-1 to 
126-MN will be described later. 

Selection combining section 128 selects M digital 
baseband signals repeatedly from among outputs of channel 

15 processing sections 126-1 to 126-MN to combine, and 
thereby forms N digital baseband signals corresponding 
to the number of reception antennas. 

Reception analog . adjusting sections 129 are 
provided corresponding to the number (N) of digital 

20 baseband signals, and compensate N digital baseband 
signals for changes in transmission characteristics 
occurring due to fluctuations in performance of N analog 
BB processing sections 52, radio circuits 53 and analog 
circuits 131 of development apparatus (reception system) 

25 50. Specific configurations of reception analog 
adjusting sections 129 will be described later. 

The digital baseband signal output from reception 
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analog adjusting section 129 is input to output interface 
section 130. In evaluating channel characteristics of 
digital BB processing sections 4 1 and 51 , digital baseband 
signal DBout is input to digital BB processing section 
5 51 of reception system 50 via switch SW4 . Meanwhile, 
in evaluating channel characteristics of digital BB 
processing sections 41 and 51 and analog BB processing 
sections 42 and 52, analog baseband signal ABout obtained 
in analog circuit 131 is input to analog BB processing 

10 section 52 of reception system 50 via switch SW3. Further, 
in evaluating channel characteristics of digital BB 
processing sections 41 and 51, analog BB processing 
sections 42 and 52, and radio circuits 43 and 53, radio 
signal RFout obtained in analog circuit 131 is input to 

15 radio circuit 53 of reception system 50. 

(2-2) Configurations of the reference channel path 
control section and channel processing sections 

Described below are configuration examples of 
reference channel path control section 127 and channel 

20 processing sections 126-1 to 126-MN. Herein, two 
configuration examples will be described in the case of 
using the expanded Sasaoka method and the case of using 
eigenvalue transformation method. It is assumed using 
parameters P10 to P20 and P30 as shown in FIG. 23 in 

25 following descriptions. 

(2-2-1) Case of using the expanded Sasaoka method 
FIG. 24 illustrates a configuration of reference 



46 



,„ Re ference channel path 
channel path control section 127. Refer 

control section 12 , is comprised of reference channel 
mod el forming section 140 and Instantaneous variatron 
initial value generating section 141. 

Beference channel model forming section 140 has 
standard ..del generatrn, section 142 to manually set 
complex i^lse response information U... set using 

gtatistical model generating 

control apparatus 121), 

section 143 that periodically updates and sets complex 
impulse response using random numhers, and actual run 
model generating section 144 that reads complex impulse 
response information obtained from Kay-Trace simulation, 
sctual running experiments and so on to seguentially 
upd ate and set, and selects in selection section 14 a 
5 channel model of a single channel generated in either 
■ o£ t he model generating sections 142 to 144 to output. 

By this means, reference channel model forming 
section 140 forms the complex impulse response 
formation .including the numher of paths and delay and 
,0 complex gain of each path, of the reference channe 1 the t 
va rie, at intervals of several docenmeters. 
each of the model generating sections 142 to 144 is of 

anH descriptions thereof are 
well-known techniques, and descrip 

omitted herein. 

instantaneous variation initial value generate 
section 141 generates an instantaneous variation mitral 
valu e of each path of the reference channel to he a random 
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value using random numbers. In addition, control 
apparatus 121 inputs to reference channel model 
generating section 140 parameter P10 (model type 
instruction to instruct a run model to select, traveling 
5 speed and direction, arrangements and di re c t i ona 1 i t i e s 
of transmission and reception antennas, and ON /OFF 
instruction of phase variation). Further, control 
apparatus 121 inputs parameter Pll (the number of paths 
and delay and complex gain of each path) to standard model 
10 generating section 142. Furthermore, control apparatus 
121 inputs parameter P12 (Ray-Trace/actual running 
experiment data) to actual run model generating section 
144 . 

As a reference channel model, selecting section 145 
15 output s parameter PI 4 ( car rier frequency , t rave 1 ing speed 
and direction, arrangements and direct ionalit ies of 
transmission and reception antennas, and ON/OFF 
instruction of phase variation), and parameter P15 (the 
number of path divisions (in compression), the number 
20 of paths of the reference channel, and delay, short-term 
variation complex gain, the angle of arrival and the 
line-of -s ight angle of each path of the reference 
channel ) . 

FIG. 25 illustrates a configuration of each of 
25 channel processing sections 126-1 to 126-MN. The 
configuration of each of channel processing sections 
126-1 to 126-MN is the same, and described below is the 
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configuration of channel processing section 126-1. In 
channel processing section 126-1, parameters P14 and P15 
are input to as s i gned- channe 1 short-term complex impulse 
response generating section 150. 
5 Assigned-channel short-term complex impulse 

response generating section 150 calculates a difference 
in distance between the reference channel and the assigned 
channel from the arrangements of transmission and 
reception antennas, calculates the complex gain of the 

10 short-term variation of each path of the assigned channel 
based on the difference in distance, and outputs the gain 
to data interpolation section 151 as parameter P18, while 
outputting the number of paths of the assigned channel, 
the delay, the angle of arrival and the line-of -sight 

15 angle of each path to correlated gaussian noise generating 
section 152 as parameter P20. 

In other words, assigned-channel short-term complex 
impulse response generating section 150 assumes that 
gains in long-term variations and short-term variations 

20 of each path included in short-term complex impulse 
response are equal in an area where transmission and 
reception antennas are installed, thereby further assumes 
that the assigned-channel and reference channel have the 
same number of paths, and that only the delay and phase 

25 of each path differ by a difference in distance obtained 
from transmission and reception points of the reference 
channel and assigned channel, a positional relationship 
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between transmission and reception antennas of the 
assigned channel, and the radiation direction and the 
direction of arrival of each path, and thus, generates 
complex impulse response of the assigned channel. 
5 More specifically, path forming section 190 

(FIG. 30) described later generates the delay, and 
assigned-channel short-term complex impulse response 
generating section 150 generates the complex gain with 
levels of the I component and Q component controlled 
10 corresponding to variations in phase . 

Data interpolation section 151 performs data 
interpolation on the complex impulse response, thereby 
performing upconverting , and outputs the resultant to 
short-term variation adding section 155 in fading adding 
15 section 154. In this way, channel processing section 
126-1 is provided with data interpolation sect ion 151, 
and therefore, even when the processing operation prior 
to data interpolation section 151 is slow to some extent, 
is capable of adding fine variations corresponding to 
20 sampling frequency f s of a baseband signal. The same 
relationship exists between data interpolation section 
153 and correlated gaussian noise generating section 152 . 

Correlated gaussian noise generating section 152 
receives parameters P14, P15 and P20, and generates the 
25 correlated gaussian noise of each path of the assigned 
channel. In other words, in channel simulator 120 of 
this Embodiment, correlated gaussian noise generating 
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sections 152 respectively in channel processing sections 
126-1 to 126-MN form correlated instantaneous variations 
of MN channels having correlation between channels, or 
between channels and paths. 
5 Correlated instantaneous variation P16 (including 

information of the number of paths and delay of each path, 
as well as the complex gain in instantaneous variation 
of each path) is interpolated in data interpolation 
section 153, and output to correlated instantaneous 
10 variation adding section 156. In addition, the 

information of the number of paths and delay of each path 
is used as information to form multipath having delays 
corresponding to the antenna arrangement as described 
later . 

15 FIG. 26 illustrates a configuration of correlated 

gaussian noise generating section 152. Correlated 
gaussian noise generating section 152 generates multitone 
having an initial phase of an initial value of the 
instantaneous value of each path of the reference channel, 

20 as the gaussian noise with a frequency band corresponding 
to the angle of arrival and 1 ine - o f - s i gh t angle of each 
path of the reference channel, performs weighting on the 
multitone with a Doppler filter and correlated filter 
characteristics using the antenna arrangement 

2 5 information as a parameter, and thus , forms the correlated 
instantaneous variation correlated with the 
instantaneous variation of the reference channel. In 



other words, the Sasaoka method as described above is 
applied . 

More specifically, multitone generating section 161 
generates multitone having as an initial phase an initial 
5 value of the instantaneous variation of each path of the 
assigned channel generated in instantaneous variation 
initial value generating section 160. The multitone is 
limited to a predetermined band within Doppler frequency 
f D by Doppler filter 162, and output to filter 165A with 

10 filter characteristics of equation (12) . 

Meanwhile, multitone generating section 163 
generates mul t itone having an initial phase corresponding 
to an initial value of the instantaneous variation of 
each path of the reference channel generated in 

15 instantaneous variation initial value generating section 
141 (FIG. 24). The multitone is limited to a 

predetermined band within Doppler frequency f D by Doppler 
filter 164, and output to filter 165B with filter 
characteristics of equation (11). 

20 At this point, Doppler filters 162 and 164 have the 

carrier frequency and traveling speed and direction input 
thereto, corresponding to which, characteristics of 
Doppler filters 162 and 164 are determined. Further, 
correlation filter section 165 has the carrier frequency, 

25 traveling speed and direction, arrangements and 
directionalities of transmission and reception antennas, 
the angle of arrival and line-of -sight angle of each path 
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input thereto, corresponding to which, characteristics 
of each of the filters 165A and 165B are determined. 

An output from correlation filter section 165 is 
added in adder 166, and input to phase variation ON /OFF 
5 section 167. Phase variation ON/OFF section 167 controls 
ON /OFF of phase variation of the correlated gaussian noise 
corresponding to an instruction on phase variation ON/OFF 
from control apparatus 121. More specifically, when 
being instructed to control the phase variation ON, the 

10 section 167 outputs the correlated gaussian noise input 
from adder 166 without any processing. 

Meanwhile, when being instructed to control the 
phase variation OFF, the section 167 obtains envelop 
amplitude V~ (I 2 + Q 2 ) of variation values of the I channel 

15 and Q channel of the correlated gaussian noise, and outputs 
the obtained envelop amplitude of variation values as 
signals of the I channel and Q channel. In other words, 
the correlated gaussian noise with the I channel and Q 
channel of the same level is formed as an instantaneous 

20 variation value, and it is thereby intended that a 
subsequent section, correlated instantaneous variation 
adding section 156, only adds a level variation without 
adding a phase variation. The reason will be described 
later. 

25 An output of phase variation ON/OFF section 167 is 

output to correlated instantaneous variation adding 
section 156 via delay section 168, as the instantaneous 
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variation of the assigned channel. 

Thus, correlated gaussian noise generating section 
152 provided for each channel obtains a correlated 
instantaneous variation correlated with the 
5 instantaneous variation of the reference channel, and 
it is thereby possible to form correlated instantaneous 
variations of MxN channels correlated with the reference 
channel from the information of each path of the reference 
channel. As a result, it is possible to simulate 

10 instantaneous variations of MxN channels with accuracy 
and ease, as compared with the case of independently 
setting the instantaneous variations of MxN channels. 

In addition, the case is described herein of 
obtaining correlated instantaneous variations 

15 correlated with the reference channel using multitone. 
Further, it may be possible that multitone generating 
sections 161 and 163 s imp ly generate white gaussian noises, 
Doppler filters 162 and 164 are set for filter 
characteristics to only pass a band in consideration of 

20 the direction of arrival of a path, and thereby correlated 
instantaneous variations of MxN channels are obtained. 

In other words, without being limited to the case 
of generating multitone and generating correlated 
instantaneous variations correlated with instantaneous 

25 variation of the reference channel, it may be possible 
to form correlated instantaneous variations of MxN 
channels by generating band-limited gaussian noises of 
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the reference channel and of another channel, performing 
weighted addition on the two band-limited gaussian noises 
with correlated filter characteristics having at least 
antenna arrangement information as a parameter, thereby 
5 forming a correlated instantaneous variation correlated 
with the instantaneous variation of the reference channel, 
and performing the aforementioned processing 
corresponding to MxN channels. 

(2-2-2) Case of using the eigenvalue 

10 transformation method 

With reference to FIGs . 27 , 28 and29, des cr ibed below 
are configurations of the reference channel path control 
section, channel processing section and correlated 
gaussian noise generating section in the case of using 

15 the eigenvalue transformation method. 

FIG. 27 illustrates a configuration of reference 
channel path control section 170 (that corresponds to 
reference channel path control section 127 in FIG. 22) 
in the case of using the eigenvalue trans formation method, 

20 with the same sections as in FIG. 24 assigned the same 
reference numerals. Unit eigenvector calculating 

section 171 as transformation matrix calculating means 
receives the information of arrangements and 
directionalities of transmission and reception antennas 

25 and the information of the angle of arrival and 
line-of -sight angle of each path of the reference channel 
among parameters P14 and P15 output from reference channel 
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model forming section 140. 

Unit eigenvector calculating section 171 first 
obtains a correlation matrix from the positional 
relationship of transmission and reception antennas, the 
5 radiation direction and angle of arrival of a signal of 
the reference channel and theoretical correlation values 
Rayleigh fading (using the spatio correlation function 
of eq. (1) in only obtaining the correlation between 
channels, while using the spatio-temporal correlation 

10 function of eq. (2) in obtaining the correlation between 
channels and paths) . The matrix has MxN rows and MxN 
columns when obtaining a correlation matrix between 
channels, while having MxNx (the number of paths) rows 
and MxNx (the number of paths) columns in obtaining a 

15 correlation matrix between channels and paths. 

As described in item (1-3-1), unit eigenvector 
calculating section 171 calculates a unit eigenvector 
(actually, which is conjugate complex transpose of the 
unit eigenvector) based on equations (3), (4) and (5), 

20 and then, outputs the eigenvector to correlated gaussian 
noise generating section 173 as a transformation matrix 
to calculate mutually correlated signal vectors from 
non- cor related signal vectors. Actually, unit 

eigenvector calculating section 171 generates an initial 

25 value of the instantaneous variation of each path of each 
channel as well as the unit eigenvector, and outputs these 
values as parameter P30 to correlated gaussian noise 
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generating section 173 in channel processing section 172 
as shown in FIG. 28. 

A configuration of channel processing section 172 
in FIG. 28 is the same as the configuration as shown in 
5 FIG. 25 except correlated gaussian noise generating 
section 173 having a different configuration, and only 
the configuration of correlated gaussian noise generating 
section 173 will be describedherein. FIG. 29 illustrates 
the configuration of correlated gaussian noise gene rating 

10 section 173 . 

Correlated gaussian noise generating section 173 
generates the number (MxNx (the number of paths) of 
instantaneous variations mutually independent between 
channels and paths in Doppler filter 180. More 

15 specifically, an initial value of the instantaneous 
variation of each path of channel 1-1 is input to 
band-limited white gaussian noise generating section 
(LWGN) 181-1, an initial value of the instantaneous 
variation of each path of channel 1-2 is input to 

20 band-limited white gaussian noise generating section 
181-2, similar processing is repeated subsequently, and 
finally, an initial value of the instantaneous variation 
of each path of channel M-N is input to band-limited white 
gaussian noise generating section 181-MN, whereby 

25 band-limited white gaussian noise generating sections 
181-1 to 181-MN generate mutually independent 
band-limited white gaussian noises. The mutually 
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independent band- 1 imi ted white gaussian noises are 
limited to bands within Doppler frequency f D in respective 
Doppler filters 182-1 to 182-MN, and then, output to 
weighted addition section 183. 
5 Weighted addition section 183 as matrix operation 

means performs the matrix operation processing using an 
eigenvector of the assigned channel on the MxNx (the number 
of paths) instantaneous variations mutually independent 
between channels and paths obtained in Doppler filer 180, 

10 and thereby obtains correlated instantaneous variations 
mutually correlated between paths. In addition, the 
correlated instantaneous variations have the correlation 
also between channels. 

The correlated instantaneous variations output from 

15 weighted addition section 183 are output to correlated 
instantaneous variation adding section 156 (FIG. 28) via 
phase ON /OFF section 184, as respective instantaneous 
variations of paths of the a s s i gned- channe 1 . 

Thus, MxNx (the number of paths) instantaneous 

20 variations are generated which are mutually independent 
between channels and paths, a correlation matrix 
( (MxNx (the number of paths) x ((MxNx (the number of paths ) ) 
is obtained from input data or experiment data, a 
difference in propagation path distance of each path is 

25 obtained from the positional relationship of antennas, 
and theoretical temporal-spatio correlation values of 
Rayleigh fading, a transformation matrix is obtained 



' • 58 

based on the correlation matrix to calculate mutually 
correlated signal vectors from mutually non-correlated 
signal vectors, and MxNx (the number of paths) 
instantaneous variations are subjected to the matrix 
5 operation processing using the transformation matrix, 
thereby obtaining MxNx (the number of paths) correlated 
instantaneous variations mutually correlated between 
paths. Therefore, it is possible to obtain correlated 
instantaneous variations mutually correlated between 

10 channels and paths, and it is thereby possible to perform 
channel simulation on MxN channels with multipath 
accurately and readily. 

Similarly, MxNx ( the numbe r o f paths ) instantaneous 
variations are generated which are mutually independent 

15 between channels, a correlation matrix (MxN x MxN) is 
obtained from input data or experiment data, a difference 
in propagation path distance of each path is obtained 
from the positional relationship of antennas, and 
theoretical spatio correlation values of Rayleigh fading, 

20 a transformation matrix is obtained based on the 
correlation matrix to calculate mutually correlated 
signal vectors from mutually non-correlated signal 
vectors, and the plurality of instantaneous variations 
are subjected to the matrix operation processing using 

25 the transformation matrix the number of times 
corresponding to the number of path, thereby obtaining 
MxN-channel correlated instantaneous variations 
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mutually correlated between channels. By this means, 
it is possible to obtain correlated instantaneous 
variations of MxN channels having correlation between 
all the channels, as well as correlation between the 
5 reference channel and each channel. As a result, it is 
possible to form instantaneous variations closer to 
instantaneous variations occurring in actual MxN channel 
transmission . 

In addition, while the case is descried herein of 

10 forming correlated instantaneous variations using the 
eigenvalue transformation method, the similar 
configuration enables formation of correlated 
instantaneous variations using the Cholesky 
factorization method as described above. 

15 As explained briefly, instead of calculating the 

eigenvector, as described in item (1-3-2), unit 
eigenvector calculating section 171 in FIG. 27 performs 
Cholesky factorization on the path correlation matrix 
to obtain a lower triangular matrix, calculates a 

20 conjugate complex transpose matrix of the matrix, based 
on equations (6) and (7), and outputs the resultant to 
correlated gaussian noise generating section 173 in 
channel processing section 172. 

In correlated gaussian noise generating section 173, 

25 the transformation matrix obtained by Cholesky 
factorization is input to weighted addition section 183 
as matrix operation means , weighted addition is performed 
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using the transformation matrix, and correlated 
instantaneous variations are thereby obtained. By this 
means, weighted addition section 183 performs the 
operation using the transformation matrix with half 
5 elements thereof of zero, and thus is capable of obtaining 
correlated instantaneous variations with a small amount 
of calculations . 

(2-3) Configuration of the fading adding section 
FIG. 30 illustrates a configuration of the fading 

10 adding section provided in each of channel processing 
sections 126-1 to 126-MN. Fading adding section 154 
receives a digital baseband signal output from switch 
125 (FIG. 22) in path forming section 190 comprised of 
shift register 191 and selector 192, and forms a signal 

15 of each path using path forming section 190. More 
specifically, shift register 191 shifts the input digital 
baseband signal by time obtained by dividing a maximum 
delay time of the path by a sampling cycle in analog BB 
processing section 42 (FIG. 21) . 

20 Selector 192 selects and outputs signals 

corresponding to the number of paths from among respective 
signals output from shift stages of shift register 191. 
At this point, path forming section 190 receives as its 
inputs the number of paths instructed from control 

25 apparatus 121 and parameter Pll indicative of delay time 
corresponding to the arrangements of transmission and 
reception antennas on a signal of each channel, and based 
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on Pll, shift register 191 and selector 192 operate. In 
this way, selector 192 in path forming section 190 outputs 
a signal of each path provided with a path delay on the 
assigned channel corresponding to the arrangements of 
5 transmission and reception antennas. 

The signal corresponding to each path is output to 
respective one of complex multipliers Al to Ak in 
correlated instantaneous variation adding section 156. 
Further, each of complex multipliers Al to Ak is supplied 
10 with correlated gaussian noise P17 output from data 
interpolation section 153. By this means, each of 
complex multipliers Al to Ak outputs a signal of each 
path provided with the correlated instantaneous 
variation . 

15 The signal of each path provided with the correlated 

instantaneous variation is output to respective one of 
a plurality of complex multipliers , BltoBk, constituting 
short -term variation adding section 155 . Each of complex 
multipliers Bl to Bk is supplied with complex gain P19 

20 of the short-term variation of each path output from data 
interpolation section 151. Short-term variation adding 
section 155 thereby outputs the signal of each path 
convoluted with complex impulse response. The signal 
of each path is all added in adders CI, C2 , , thereby 

25 forming a multipath signal in which channel variations 
are reflected. 

The multipath signal is supplied to adder C3 . Adder 
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C3 is supplied with the white gaussian noise which is 
generated in white gaussian noise generating section 
(WGN) 21 and amplified in amplifier 22 to noise level 
S4 designated from control apparatus 30. Thus, adder 
5 C3 adds receiver noise to the multipath signal- 
Fading adding section 154 further has automatic gain 
control section 193. In automatic gain control section 

193, AGC control section 195 receives a target level from 
control apparatus 121, and thereby sets as an 

10 amplification value of amplifier 194 a difference value 
between the target level and an output signal of amplifier 

194. As a result, automatic gain control section 193 
performs simplified digital gain control processing, and 
is capable of making the multipath signal a constant signal 

15 with the target level. 

The reason why it is necessary to perform gain 
control on a multipath signal is because the multipath 
signal resulting from addition by adder CI is obtained 
by adding respective signals of paths each provided with 

20 the level variation independently, and it is assumed that 
the level variation occurs in the digital baseband signal . 
In consideration thereof, gain control section 193 
executes the simplified digital gain control processing 
to make the level of a multipath signal constant, and 

25 it is thereby possible to prevent a bit from being lost 
in AD conversion in reception system 50 of the development 
apparatus even when radio circuit 53 (FIG. 21) is not 
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completed and the AGC processing cannot be carried out. 
As a result, it is possible to evaluate channel 
characteristics on a multipath channel with excellence 
based on the digital baseband signal of digital BB 
5 processing section 41. 

Further, in channel simulator 120 of this 
Embodiment, a digital baseband signal is input from 
digital BB processing section 41 of the transmission 
system, the signal is provided with channel variations, 

10 and the resultant signal is output to digital BB processing 
section 51 of reception system 50. In evaluating 
propagation characteristics of digital BB processing 
sections 41 and 51, phase variation ON /OFF sections 167 
(FIG. 26) and 184 (FIG. 29) are set for OFF, and the 

15 correlated instantaneous variation with the I channel 
and Q channel of the same level is input to correlated 
instantaneous variation adding section 156. This is not 
shown in the figure, but is the same as in the short-term 
variation supplied to short-term variation adding section 

20 153. 

By this means, even when AFC of radio circuit 53 
is not available, it is possible to evaluate performance 
of digital BB processing sections 4 1 and 51 . In addition, 
in a stage where radio circuits 43 and 53 are connected, 
25 since the AGC function and AFC function due to radio 
circuit 53 serve, each of complex multipliers Al to Ak 
and Bl to Bk may multiply a complex gain of the short-term 
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variation such that envelop amplitude is different 
between the I channel and Q channel to provide a digital 
baseband signal with a phase variation. 

(2-4) Configurations of the transmission analog 
5 adjusting section and reception analog adjusting section 
Described below are configurations of transmission 
analog adjusting section 124 and reception analog 
adjusting section 129. Transmission analog adjusting 
section 124 and reception analog adjusting section 129 
10 are to simulate fluctuations in a signal of each channel 
occurring due to fluctuations in performance of analog 
circuitry corresponding to the each channel among MxN 
channels . 

In other words, attention was directed toward the 
15 fact that development apparatuses 40 and 50 targeted for 
simulation have M analog circuits on the transmission 
side and N analog circuits on the reception side, and 
that fluctuations between the MxN analog circuits affect 
signals on channels, and transmission analog adjusting 
20 section 124 and reception analog adjusting section 129 
simulate the fluctuations between the channels to add 
to a digital baseband signal as appropriate. It is 
thereby possible to simulate propagation variations in 
MxN-channel t r ansmi s s i on mor e closely approximating real 
25 variations . 

As a result, in channel simulator 200, even if 
development of radio circuit 43 of transmission system 



65 

40 and radio circuit 53 of reception system 50 has not 
been completed, it is possible to add signal deterioration 
assumed to occur in radio circuits 43 and 53 to the digital 
baseband signal to evaluate characteristics of digital 
5 BB processing sections 41 and 51. 

As a result, it is possible to evaluate channel 
characteristics of digital BB processing sections 41 and 
51 including appropriateness between digital BB 
processing sections 41 and 51 and radio circuits 43 and 
10 53. Further, it is possible to beforehand measure an 
extent of deterioration to occur in radio circuits 43 
and 53 subsequently developed with the performance of 
digital BB processing sections 41 and 51 meeting a desired 
value . 

15 Configurations of transmission analog adjusting 

section 124 and reception analog adjusting section 129 
will specifically be describedbelow. As shown in FIG. 31, 
in transmission analog adjusting section 124, a baseband 
signal from switch 125 (FIG. 22) is input to gain unbalance 

20 generating section 210. Gain unbalance generating 
section 210 amplifies the I channel signal and Q channel 
signal of the digital baseband signal independently, and 
thereby generates a gain difference. DC offset adding 
section 211 increases or decreases each of the I channel 

25 signal and Q channel signal by a constant value, and 
thereby adds the DC offset. 

Frequency offset • phase offset adding section 212 



66 

adds the frequency offset and phase offset assumed to 
occur in radio circuit 43 and analog circuit 123 (FIG. 22) 
to the I channel signal and Q channel signal. Actually, 
frequency offset • phase offset adding section 212 is 
5 comprised of a complex multiplier which multiplies each 
channel signal by variation amount cos 9 1 or sin 9 2 
respectively corresponding to instantaneous phase 0 1 or 
9 2 . 

In other words, the I channel signal is multiplied 

10 by variation amount cos 9 1, while the Q channel signal 
is multiplied by variation amount sin 9 2. Herein, 
instantaneous phases 9 1 and 9 2 made constant means only 
adding the phase offset, while instantaneous phases 9 1 
and 9 2 varying with time means adding the frequency offset 

15 in addition to the phase offset. 

In obtaining instantaneous phases 9 1 and 9 2, 
transmission analog adjusting section 124 calculates a 
phase rotation amount per sample from frequency offset 
set value S20E in phase increment calculating section 

20 2 15 to output to mod 2 n calculating circuits 217 and 219. 
At this point, to add deterioration in orthogonality 
between the I channel signal and Q channel signal, adder 
218 adds orthogonality deterioration amount S20F to a 
phase rotation amount of the Q channel signal. 

25 A phase of the last sample is input to adder 216. 

The phase of the last sample is calculated by Z-l 
calculating circuit 222 performing computation based on 
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an initial phase (i.e. phase offset) S20D and the phase 
of the last sample. Adder 216 adds the phase rotation 
amount corresponding to one sample calculated in phase 
increment calculating circuit 215 to the phase of the 
5 last sample, and thereby obtains a phase rotation amount 
of a current sample. 

Thus, by repeating a processing loop of adder 216, 
mod 2 n calculating circuit 217 and Z-l calculating 
circuit 222, I channel instantaneous phase 0 1 provided 

10 with the phase offset and frequency offset is calculated 
for each sample, and Q channel instantaneous phase Q 2 
is calculated by adding the deterioration amount of 
orthogonality to I channel instantaneous phase 8 1. 

Then, frequency offset • phase offset adding section 

15 212 adds variation amount cos 9 1 and variation amount 
sin 6 2 respectively to the I channel and Q channel of the 
digital baseband signal, and thereby adds the frequency 
offset and phase offset on each channel of the digital 
baseband signal assumed to occur in radio circuit 43 of 

20 transmission system 40 and analog circuit 123. Delay 
adjusting section 213 adds a circuit delay amount assumed 
to occur in radio circuit 43 and the analog circuit. 

Pseudo power amplifier (PA) section 214 is to 
generate simulated non-linear distortion assumed to occur 

25 in an amplifying section of radio circuit 43, and for 
example, is configured as shown in FIG. 32. Pseudo PA 
section 214 calculates V" (I 2 + Q 2 ) in envelop amplitude 
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calculating section 230, and thereby calculates envelop 
amplitude X of the digital baseband signal to output to 
averaging circuit 231 and distortion calculating section 

232 . 

5 Averaging circuit 231 averages the envelop 

amplitude for a time corresponding to forgetting factor 
(i.e. level calculation time constant ) S20H set by control 
apparatus 121, and outputs obtained average value Pave 
to saturation level computation circuit 233. Saturation 
10 level- computation circuit 233 obtains saturation level 
Asat in the following equation, as suming the average value 
of envelop amplitude as Pave and backoff of the power 
amplifier set in control apparatus 121 as IBO. 

1BO 

A sat =P ave xi0' 20 (13) 

15 

Distortion computation section 232 calculates a 
control value of amplifier 2 34 in the following equation, 
using the envelop amplitude value X obtained in envelop 
amplitude calculating circuit 230 and saturation level 
20 Asat obtained in saturation level computation circuit 

233 . 

Control 1 



value J_ 



1*1 



A sat 



To (14) 



In this way, pseudo power amplifier (PA) section 
214 is capable of adding simulated non-linear distortion 
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assumed to occur in the amplifying section of radio circuit 
43 to the digital baseband signal. 

Reception analog adjusting section 129 is 
configured as shown in FIG. 33. In reception analog 
5 adjusting section 129, the digital baseband signal output 
from selection combining section 128 (FIG. 22) is input 
to frequency offset ■ phase offset adding section 251. 

Frequency offset * phase offset adding section 251 
performs the same processing as in frequency offset ■ 

10 phase offset adding section 212 of transmission analog 
adjusting section 124. In other words, the section 251 
adds the frequency offset and phase offset assumed to 
occur in radio circuit 53 of reception circuit 50 and 
analog circuit 131 (FIG. 22) t o I and Q channe 1 s . Actually, 

15 frequency offset • phase offset adding section 251 is 
comprised of a complex multiplier which multiplies each 
channel signal by variation amount cos 0 1' or sin 0 2' 
respectively corresponding to instantaneous phase 0 1' 
or 0 2'. In other words, the I channel signal is 

20 multiplied by variation amount cos 0 1', while the Q 
channel signal is multiplied by variation amount sin 0 2'. 

In obtaining instantaneous phases 0 1' and 0 2' , 
reception analog adjusting section 12 9 calculates a phase 
rotation amount per sample from frequency offset set value 

25 S22B in phase increment calculating section 252 to output 
to mod 2 7i calculating circuits 254 and 256. At this 
point, to add deterioration in orthogonality between the 
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I channel signal and Q channel signal, adder 255 adds 
orthogonality deterioration amount S22C to a phase 
rotation amount of the Q channel signal. 

A phase of the last sample is input to adder 253. 
5 The phase of the last sample is calculated by Z-l 
calculating circuit 259 performing computation based on 
an initial phase (i.e. phase offset) S22A and the phase 
of the last sample. Adder 253 adds the phase rotation 
amount corresponding to one sample calculated in phase 

10 increment calculating circuit 252 to the phase of the 
last sample, and thereby obtains a phase rotation amount 
of a current sample. 

Then, frequency offset • phase offset adding section 
251 adds variation amount cos d 1' and variation amount 

15 sin 6 2' respectively to the I channel and Q channel of 
the digital baseband signal, and thereby adds the 
frequency offset and phase offset on each channel of the 
digital baseband signal assumed to occur in radio circuit 
53 of reception system 50 and analog circuit 131. 

20 Gain unbalance generating section 261 amplifies the 

I channel signal and Q channel signal of the digital 
baseband signal independently, and thereby generates a 
gain difference. DC offset adding section 2 62 increases 
or decreases each of the I and Q channels by a constant 

25 value, and thereby adds the DC offset. Delay adjusting 
section 263 adds a circuit delay amount assumed to occur 
in radio circuit 53 and analog circuit 131. 
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A user is capable of selecting arbitrarily via 
control apparatus 121 various set values S20 (S20A to 
S20I) and S22 (S22A to S22H) of transmission analog 
adjusting section 124 and reception analog adjusting 
5 section 12 9. 

Thus, it is possible to simulate gain unbalance, 
DC offset, frequency offset, phase offset, circuit delay , 
non-linear distortion in amplification, and so on assumed 
to occur in radio circuits 43 and 53 and analog circuits 
10 123 and 131, prior to completion of radio circuit 43 of 
transmission system 40 and radio circuit 53 of reception 
system 50, i.e. at the stage where digital BB processing 
sections 41 and 51 have only been completed . Itisthereby 
possible to evaluate characteristics of digital BB 
15 processing sections 41 and 51 in combining digital BB 
processing sections 41 and 51 under development and radio 
circuits 43 and 53 with various characteristics. 
(3) Effects of the Embodiment 

Thus, according to the aforementioned constitution, 
20 by providing switch 125 which makes N copies of each of 
M signals obtained in transmission system 40, and thereby 
forms MxN channel signals, channel processing sections 
126-1 to 126-MN which respectively add correlated 
instantaneous variations and short-term variations 
25 corresponding to arrangements of 'transmission and 
reception antennas to the MxN channel signals, and 
selection combining section 128 which selectively 
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combines M signals repeatedly among the MxN signals 
provided with channel variations to form N signals, it 
becomes possible to simulate channel variations actually 
occurring in a mu 1 1 i - ant enna apparatus, and it is thereby 
5 possible to simulate channel characteristics in the 
multi-antenna apparatus with accuracy and ease. 

The present invention is not limited to the 
aforementioned Embodiment, and is capable of being 
carried into practice with various modifications thereof . 

10 A channel simulation method according to an aspect 

of the invention includes a channel variation forming 
step of forming a channel variation on each of MxN channels 
using information of arrangements of transmission and 
reception antennas, and a channel variation adding step 

15 of adding channel variations corresponding to the MxN 
channels to respective signals of the MxN channels. 

According to this method, the channel variations 
of all the MxN channels are formed from the information 
of arrangements of transmission and reception antennas, 

20 and it is thereby possible to form channel variations 
in MxN-channel transmission formed by the multi-antenna 
apparatus with accuracy and ease. Further, in the case 
of collecting channel data in running experiments to 
reproduce, data corresponding to a single channel is 

25 collected by a data collection apparatus having a 
transmission antenna and a reception antenna, and using 
the data as reference channel data, channel variations 
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in MxN-channel transmission can be formed accurately and 
readily from the reference channel data and relative 
arrangements of transmission and reception antennas of 
a development apparatus. Thus, the memory for data 
5 storage can be saved largely, the number of times of 
running experiment thereby decreases drastically, and 
it is possible to improve efficiency in development. 

In a channel simulation method of another aspect 
of the invention, in the channel variation forming step, 

10 a delay and a phase variation on each channel due to the 
arrangements of antennas are obtained using the 
information of arrangements of transmission and reception 
antennas, and channel variations are formed such that 
the delay and the phase variation vary with the channels. 

15 According to this method, channels variations of 

the MxN channels are formed by only varying the delay 
and the phase variation on each of the channels caused 
by the arrangements of the antennas, and it is thereby 
possible to form the channel variations of the MxN channels 

20 with ease . 

In a channel simulation method of still another 
aspect of the invention, in the channel variation forming 
step, in forming a short-term variation on each channel 
as a channel variation, short-term variations 

25 corresponding to MxN channels are formed by obtaining 
a difference in path distance between each path of a 
reference channel beforehand set or prepared and 
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pertinent each path of each channel using information 
of a positional relationship between transmission and 
reception antennas on each channel and information of 
a radiation direction and a direction of arrival on each 
5 path, and for a signal of pertinent each path of each 
channel, generating a short-term variation such that a 
phase difference occurs with respect to a short-term 
variation of each path of the reference channel by the 
difference in path distance. 

10 According to this method, by regarding the number 

of paths on each channel and gains on paths as being equal 
since a distance between transmission antennas and a 
distance between reception antennas are sufficiently 
smaller than a short-term variation period, the 

15 short-term variation is formed such that a phase 
difference occurs with respect to the short-term 
variation of each path of the reference channel by the 
difference in path distance, whereby it is possible to 
form short-term variations of all MxN channels from a 

20 channel model of the reference channel, and thus, 
preparing beforehand only the channel model of the 
reference channel enables the short-term variations of 
the MxN channels to be formed with ease and accuracy. 

In a channel simulation method of a further aspect 

25 of the invention, in the channel variation forming step, 
in forming an instantaneous variation on each channel 
as a channel variation, correlated instantaneous 
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variations corresponding to MxN channels are formed by 
repeating processing, the number of times corresponding 
to MxN channels, for generating respective band-limited 
gaussian noises corresponding to the reference channel 
5 and another channel, subjectingtwo band- limited gaussian 
noises to weighted addition with correlated filter 
characteristics using at least the information of 
arrangements of antennas as a parameter, and thereby 
forming a correlated instantaneous variation correlated 
10 with the instantaneous variation on the reference 
channe 1 . 

According to this method, it is possible to form 
MxN-channel correlated instantaneous variations 
correlated with the reference channel from the 

15 information of each path of the reference channel, and 
it is thus possible to form instantaneous variations of 
MxN channels with accuracy and ease, as compared to the 
case of independently setting instantaneous variations 
of MxN channels. In other words, this method is to expand 

20 the method by Sasaoka conventionally proposed as a method 
of generating two-channel correlated instantaneous 
variations so as to generate MxN-channel correlated 
instantaneous variations . 

In a channel simulation method of a still further 

25 aspect of the invention, the channel variation forming 
step includes the steps of generating MxNx (the number 
of paths) instantaneous variations mutually independent 
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between channels, obtaining an (MNxMN) co r r e 1 a t i on ma t r ix 
from a difference in propagation path distance of each 
path obtained from input data or experiment data and the 
positional relationship of antennas, and theoretical 
5 spatio correlation values of Rayleigh fading, obtaining 
based on the correlation matrix a transformation matrix 
to calculate mutually correlated signal vectors from 
signal vectors that are not correlated with one another, 
and obtaining MxNx(the number of paths) correlated 

10 instantaneous variations correlated between channels, 
by repeating, the number of times corresponding to the 
number of paths, matrix operation processing using the 
transformation matrix for each instantaneous variation 
of a pertinent path of each channel. 

15 According to this method, as well as the correlation 

between the reference channel and each channel, it is 
possible to form MxNx (the number of path) correlated 
instantaneous variations correlated between all the 
channels. As a result, it is possible to form 

20 instantaneous variations closer to instantaneous 
variations occurring in actual MxN-channel transmission . 

In a channel simulation method of a yet further 
aspect of the invention, the channel variation forming 
step includes the steps of generating MxNx (the number 

25 of paths) instantaneous variations mutually independent 
between channels and between paths , obtaining an (MNx(the 
number of pa ths ) xMNx ( the number of paths)) correlation 
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matrix from a difference in propagation path distance 
of each path obtained from input data or experiment data 
and the positional relationship of antennas, and 
theoretical temporal-spatio correlation values of 
5 Rayleigh fading, obtaining based on the correlation 
matrix a transformation matrix to calculate mutually 
correlated signal vectors from signal vectors that are 
not correlated with one another, and obtaining MxNx (the 
number of paths) correlated instantaneous variations 

10 correlated between the paths, by performing matrix 
operation processing using the transformation matrix on 
the MxNx (the number of paths) instantaneous variation. 

According to this method, correlated instantaneous 
variations are obtained that are correlated between 

15 channels and between paths, and it is thereby possible 
to also perform simulations in MxN- channe 1 multipath 
transmission accurately and readily. 

In a channel simulation method of a yet further 
aspect of the invention, in the step of obtaining a 

20 transformation matrix, the transformation matrix is 
obtained by eigenvalue transformation. 

According to this method, in obtaining correlated 
instantaneous variations from mutually independent 
instantaneous variations, it is possible to use a matrix 

25 (eigenvalue) with a small number of elements, instead 
of using a matrix with (MxN) 2 or (MxNx (the number of 
paths ) ) 2 elements , and thus, the correlated instantaneous 
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variations can be obtained with a small amount of 
calculations . 

In a channel simulation method of a yet further 
aspect of the invention, in the step of obtaining a 
5 transformation matrix, the transformation matrix is 
obtained by Cholesky factorization. 

According to this method, in obtaining correlated 
instantaneous variations from mutually independent 
instantaneous variations, a lower triangular matrix 
10 obtained by Cholesky factorization is used, instead of 
using a matrix with (MxN) 2 or (MxNx (the number of paths) ) 2 
elements, and thus, the correlated instantaneous 
variations can be obtained with a small amount of 
calculations . 

15 A channel simulator of an aspect of the invention 

is a channel simulator that simulates channel 
characteristics of a wireless apparatus using an 
MxN-channel transmission system using M transmission 
antennas andN reception antennas , and has a configuration 

20 provided with an input section which inputs M signals 
obtained by a transmission system of the wireless 
apparatus, a signal replicating section which makes N 
copies of each of the M signals, and thereby forms MxN 
channel signals, a channel processing section that adds 

25 . a channel variation to each of the MxN channel signals 
corresponding to arrangements of transmission and 
reception antennas, and a combining section that 
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selectively combines M channel signals repeatedly among 
the MxN channel signals each provided with the channel 
variation to form N signals. 

According to this configuration, it is possible to 
5 simulate channel variations actually occurring in a 
multi-antenna apparatus, and it is thus possible to 
simulate channel characteristics in the multi-antenna 
apparatus with accuracy and ease. 

In a channel simulator of another aspect of the 

10 invention, the channel processing section is provided 
with a path forming section that forms a signal of each 
path having a delay corresponding to the arrangements 
of transmission and reception antennas for a signal of 
each channel, a short-term complex impulse response 

15 generating section that forms a complex gain of a 
short-term variation to be added to each path of each 
channel, and a short-term variation adding section that 
adds the short-term variation to the signal of each path 
of the each channel, and the short-term complex impulse 

20 response generating section obtains a difference in path 
distance between each path of a reference channel and 
pertinent each path of each channel using information 
of a positional relationship between transmission and 
reception antennas on each channel and a radiation 

25 direction and a direction of arrival on each path, and 
for the signal of each path of each channel generated 
in the path forming section, generates a short-term 
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variation such that a phase difference occurs with respect 
to a short-term variation of each path of the reference 
channel beforehand set or prepared by the difference in 
path distance . 

5 According to this configuration, it is possible to 

form short-term variations of all MxN channels from a 
channel model of the reference channel, and thus, 
preparing beforehand only the channel model of the 
reference channel enables the short-term variations of 

10 the MxN channels to be formed with ease and accuracy. 

A channel simulator of still another aspect of the 
invention adopts a configuration where the channel 
processing section is provided with a path forming section 
that forms a signal of each path having a delay 

15 corresponding to the arrangements of transmission and 
reception antennas for a signal of each channel, a 
correlated gaussian noise generating section that 
generates a correlated instantaneous variation to be 
added to each path of each channel, and a correlated 

20 instantaneous variation adding section that adds the 
. correlated instantaneous variation to the signal of each 
path of each channel. 

According to this configuration, as compared with 
the case of independently setting MxN-channel 

25 instantaneous variations, it is possible to excellently 
simulate MxN-channel instantaneous variations assumed 
to actually occur in a multi-antenna apparatus to add 
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to the signal of each path provided with the delay 
corresponding to the arrangements of transmission and 
reception antennas . 

A channel simulator of a further aspect of the 
5 invention adopts a configuration where the correlated 
gaussian noise generating section forms correlated 
instantaneous variations corresponding to MxN channels 
by repeating processing, the number of times 
corresponding to MxN channels, for generating respective 

10 band-limited gaussian noises corresponding to the 
reference channel and another channel, subjecting two 
band-limited gaussian noises to weighted addition with 
correlated filter characteristics using at least the 
information of arrangements of antennas as a parameter, 

15 and thereby forming a correlated instantaneous variation 
correlated with the instantaneous variation on the 
reference channel . 

According to this configuration, it is possible to 
form MxN-channel correlated instantaneous variations 

20 correlated with the reference channel from the 
information of each path of the reference channel, and 
it is thus possible to form MxN-channel instantaneous 
variations with accuracy and ease, as compared to the 
case of independently setting MxN-channel instantaneous 

25 variations. In other words, this method is to expand 
the method by Sasaoka conventionally proposed as a method 
of generating two-channel correlated instantaneous 
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variations so as to generate MxN-channel correlated 
instantaneous variations . 

A channel simulator of a still further aspect of 
the invention adopts a configuration further provided 
5 with a transformation matrix calculating section which 
obtains a correlation matrix from a difference in 
propagation path distance of each path obtained from input 
data or experiment data and the positional relationship 
of antennas and theoretical spatio correlation values 

10 of Rayleigh fading, and then, based on the correlation 
matrix, obtains a transformation matrix to calculate 
mutually correlated signal vectors from signal vectors 
that are not correlated with one another, where the 
correlated gaussian noise generating section is provided 

15 with an instantaneous variation generating section that 
generates MxNx (the number of paths) instantaneous 
variations mutually independent between channels, and 
a matrix operation section that generates MxNx (the number 
of paths) correlated instantaneous variations correlated 

20 between channels, by repeating matrix operation 
processing using the transformation matrix on the 
instantaneous variations the number of times 
corresponding to the number of paths. 

According to this configuration, as well as the 

25 correlation between the reference channel and each 
channel, it is possible to form MxNx(the number of path) 
correlated instantaneous variations correlated between 
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all the channels. As a result, it is possible to form 
instantaneous variations closer to instantaneous 
variations occurring in actual MxN - channe 1 transmi ssion . 

A channel simulator of a yet further aspect of the 
5 invention adopts a configuration' further provided with 
a trans format ion matrix calculating section which obtains 
a correlation matrix from a difference in propagation 
path distance of each path obtained from input data or 
experiment data and the positional relationship of 

10 antennas and theoretical temporal- spa t io correlation 
values of Rayleigh fading, and then, based on the 
correlation matrix, obtains a transformation matrix to 
calculate mutually correlated signal vectors from signal 
vectors that are not correlated with one another, where 

15 the correlated gaussian noise generating section is 
provided with an instantaneous variation generating 
section that generates MxNx (the number of paths) 
instantaneous variations mutually independent between 
channels and between paths, and a matrix operation section 

20 that generates MxNx (the number of paths) correlated 
instantaneous variations correlated between the paths, 
by performing matrix operation processing using the 
transformation matrix on the instantaneous variations. 

According to this configuration, correlated 

25 instantaneous variations are obtained which are 
correlated between channels and between paths, and it 
is thereby possible to simulate characteristics in 
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MxN-channel multipath transmission accurately and 
readily . 

A channel simulator of a yet further aspect of the 
invention adopts a configuration where the transformation 
5 matrix calculating section obtains a transformation 
matrix by eigenvalue transformation. 

According to this configuration, in obtaining 
correlated instantaneous variations from mutually 
independent instantaneous variations in the matrix 

10 operation section, it is possible to use a matrix 
(eigenvalue) with a small number of elements, instead 
of using a matrix with ( MxN ) 2 or (MxNx (the number of 
paths) ) 2 elements, and thus, it is possible to reduce an 
amount of calculations in the matrix operation section. 

15 A channel simulator of a yet further aspect of the 

invention adopts a configuration where the transformation 
matrix calculating section obtains a transformation 
matrix by Cholesky factorization. 

According to this configuration, in obtaining 

20 correlated instantaneous variations from mutually 
independent instantaneous variations in the matrix 
operation section, a lower triangular matrix obtained 
by Cholesky factorization is used, instead of using a 
matrix with (MxN) 2 or (MxNx (the number of paths) ) 2 e lement s , 

25 and thus, it is possible to reduce an amount of 
calculations in the matrix operation section. 

A channel simulator of a yet further aspect of the 
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invention adopts a configuration further provided with 
an analog adjusting section which is comprised of a digital 
circuit, and simulates fluctuations in a signal of each 
channel caused by fluctuations in performance of an analog 
5 circuit corresponding to each channel among the MxN 
channels . 

According to this configuration, attention was 
directed toward the fact that the multi-antenna apparatus 
targeted for simulation has M analog circuits on the 

10 transmission side and N analog circuits on the reception 
side, and that fluctuations between the MxN analog 
circuits affect signals on channels, and the analog 
adjusting section simulates the fluctuations between the 
channels to add to a digital baseband signal as appropriate. 

15 It is thereby possible to simulate channel variations 
in MxN-channel transmission more closely approximating 
real variations. 

A channel simulator of a yet further aspect of the 
invention adopts a configuration further provided with 

20 an input interface that inputs an output signal of a 
digital baseband processing section of a transmission 
system of the wireless apparatus, a gain control section 
that performs gain control such that a signal level becomes 
almost constant of a multipath signal resulting from 

25 addition of the signal of each path provided with the 
channel variation, and an output interface that outputs 
the digital baseband signal subjected to the gain control 
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to a digital baseband processing section of the reception 
system of the wireless apparatus, where the channel 
processing section adds a channel variation component 
with an I component and a Q component equal to each other. 
5 According to this configuration, the digital 

baseband signal is directly input from the input section, 
the multipath signal provided with the channel variation 
is subjected to a level correction in the gain control 
section so as not to loose a bit in AD conversion in the 

10 reception system, and further is provided with the channel 
variation component with the I component and Q component 
equal to each another, and it is thereby possible to 
measure characteristics of the time AFC and AGC almost 
ideally operates on each path even when a radio circuit 

15 of the reception system of the development apparatus is 
not present. As a result, it is possible to evaluate 
performance of the digital baseband processing section 
only using a digital baseband signal without the AGC 
circuit and AFC circuit. It thus becomes possible to 

20 evaluate characteristics of the digital baseband 
processing section without the radio circuit, thereby 
enabling improved efficiency in development. 

As described above, according to the present 
invention, a channel variation on each of MxN channels 

25 is formed using arrangement information of reception 
antennas, and channel variations corresponding to the 
MxN channels are added to respective signals of the MxN 



87 

channels. Therefore, the channel variations of all the 
MxN channels can be formed from the information of 
arrangements of transmission and reception antennas, and 
it is thereby possible to form channel variations in 
MxN-channel transmission formed by the mul t i -ant enna 
apparatus with accuracy and ease. 

This application is based on the Japanese Patent 
Application No . 2 0 0 2 - 3 7 2 9 6 0 filed on December 24 , 2002 , 
entire content of which is expressly incorporated by 
reference herein. 

Industrial Appli cabil i ty 

The present invention is suitable for use in 
developing, for example, a cellular telephone, base 
station of the cellular telephone, and MT (Mobile 
Terminal) and AP (Access Point) of wireless LAN (Local 
Area Network) . 
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